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ABSTRACT 


We report the first detection of the ground state fine structure transition 
of Si + at a rest wavelength determined to be 34.815 ± 0.004pm. These 
observations were obtained with the facility spectrometer on NASA's Kuiper 
Airborne Observatory. A 6' NW-SE strip scan across the Orion-KL region shows 
[Sill] emission from both the extended photodissociation region surrounding ©‘Ori 
C and from the shocked gas NW of BN-KL. The inferred gas-phase silicon elemental 
abundance relative to hydrogen in the dense (~10^ cnT^) primarily neutral 
photodissociation region is approximately 2.6x10“^, a factor of 0.075 times the 
solar value and 3.4 times greater than the abundance in the moderate density (^10^ cm~3) 
cloud toward t, Oph. The silicon abundance in the shocked gas is approximately 
solar, indicating that any pre-existing grains have been destroyed in the shock 
wave or that the preshock gas carries a near solar abundance of gas phase 
silicon. The shock-excited [Sill](34.8pm) emission may arise from shocked wind 
material in the outflow around IRc2, with wind velocities £100 km/s. 

I. INTRODUCTION 

The gas-phase elemental silicon abundance in different phases of the 
interstellar medium is important for understanding interstellar chemistry and the 
formation and destruction of interstellar dust. Ultraviolet observations of 
gas-phase silicon depletions in diffuse clouds coupled with infrared observations 
of solid-state silicate absorption features suggest that a large fraction (£90$) 
of the silicon may reside in refractory silicate material in the interstellar 
medium (cf e Savage and Mathis 1979). However, theoretical models of the 
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destruction of interstellar silicate grains by supernova-driven shock waves 
indicate that interstellar silicate grains may be destroyed faster than they can 
be reinjected by red giant winds and/or supernovae ejecta (Draine and Salpeter 
1979; Dwek and Scalo 1980; Seab and Shull 1985). Consequently, in the absence of 
any other grain formation process, silicon would be expected to suffer little 
depletion in the diffuse interstellar medium. The apparent discrepancy between 
the observed depletion and the theoretical expectations could be resolved if, for 
example, silicon accreted onto grains in the denser phases of the interstellar 
medium (preferably in the form of silicates to produce the 9.7 pm feature) or if 
the previous theoretical models of shock destruction of silicate grains seriously 
overestimate the efficiency of the process. These possibilities can be 
constrained by estimates of the silicon depletion in the dense neutral gas which 
produces strong [Si II] (34. 8pm) emission. 

Silicon lies below carbon in the periodic table and has a similar electronic 
structure. The [SillDG^.Spm) line is a 2p 3 / 2 — >2p i /2 8 round state magnetic 
dipole transition analogous to the [CII](157.7pm) transition first detected by 
Russell et al. (1980). These two lines, as well as the [OI](63pm) line, are the 
dominant coolants in T <, 5000 K atomic interstellar gas which includes 
FUV-ionized (hv < 13-6 eV) metals (e.g., carbon and silicon) in the singly 
ionized state (cf Dalgarno and McCray 1972). In such a gas, with relative 
electron abundances < 10~3 f atomic hydrogen collisions dominate the excitation of 
these lines. Various line parameters are given in Table 1 ; the range in 
excitation energy and critical density illustrate that these lines probe the 
density, temperature, and elemental abundance of the gas. High temperatures (T 
100 K) and densities (>, 10 H cnT^) are prerequisites for [01] and [Sill] to 
dominate the cooling. Such high pressures are found, for example, in the 



photodissociation regions (PDRs) bordering HII regions (Tielens and Hollenbach 
1 985a, b) and behind high velocity (v > 50 km/s) shock waves which heat, 

O /sJ 

compress, dissociate, and ionize the preshock gas (cf McKee, Chernoff, and 
Hollenbach 1984). 

Previous observations of Orion-KL have demonstrated the presence of a PDR 
producing [OI](63, 145pm), [CII](157.7pm), [CI](370, 610pm), and low-J 
transitions of CO (references in Tielens and Hollenbach 1985b). Just north of 
BN-KL is a region of shocked gas which was first detected by its 2pm H 2 emission 
(Gautier et al. 1976) and has now been seen in a variety of other species 
including high-J rotational transitions of CO (Goldsmith et al. 1981; Storey et 
al. 1 981 , Watson 1982), rotational transitions , of OH (Storey, Watson, and Townes 
1981; Viscuso et al. 1985), and the ground state fine structure transition of 
[OI](63pm) (Werner et al. 1984). Here we report the first detection of the 
[SiII](34.8pm) line in a series of observations spanning Orion-KL and interpret 
the emission as arising in both the PDR and in the shocked gas. 

II. THE OBSERVATIONS 

The observations were carried out with the 91 cm telescope of the Kuiper 
Airborne Observatory on February 14, 1985 using the facility cooled grating 
spectrometer (CGS) described by Erickson et al. (1984, 1985). The beam had a 
FWHM of 46.5” with an effective area corresponding to a 48.3" diameter disc. The 
chopper throw was -4* in azimuth (approximately E-W). The resolution was 114 
km/s. The total bandpass of the six Ge:Ga photoconductor detectors was 
- 3.3 resolution elements, sufficient to include both the line and the adjacent 
continuum using a single grating position. Sequences of four 10 s integrations 



were taken with the source placed alternately in right and left beam. The 
[SiII](34.8pm) line was observed at 8 spatial positions in the Orion region which 
lie in a NW-SE line through KL and the Trapezium. In Fig. 1 we show these 
positions superimposed on the 5 GHz radio continuum map of Martin and Gull 
(1976), which shows the physical extent of the ionized gas. Absolute positions 
were determined to ±5" by offsetting from ©‘Ori C or 0 2 Ori A. The infrared 
boresight established on the ground was verified in flight to ±5", assuming that 
the 63.07pm continuum peaks at KL (R. A. (1950)=5 h 32 m 46.7 s , 6( 1 950)=-5°24 '27") 
(Werner 1982). 

The absolute flux calibration and the relative detector response were 
determined from observations of KL In the continuum adjacent to the line. The 
calibrations done at the beginning and end of the observing leg agreed to ±4$. 
This spectral region has no telluric absorption features >2$ and the instrument 
response function is flat to 10$. Hence the final reduced spectra are simply the 
signals obtained at the map positions divided by the signals obtained on the KL 
continuum and multiplied by the flux from KL as measured by Erickson et 
al. (1981) in a similar sized beam (50"). The effects of diffraction in the 
telescope were estimated to be < 3 % and were ignored. The overall uncertainty in 
the absolute calibration is probably less than ±25$. 

The spectra obtained at the ionization bar (P3)» the Trapezium (P5), and NW 
of KL (P6) are shown in Fig. 2. Since this was the first observation of the 
[SiII](34.8pm) line, we used two grating positions to confirm it at the first 
spatial position observed (P3). The 30lid curves in Fig. 2 are fits to the data 
using a flat continuum and the width, intensity, and wavelength of a superimposed 
gaussian line profile. The wavelength of the [Sill] transition is given as 
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34. 804pm by Morton and Smith (1973) and as 34.814pm by Moorwood et al. (1978). 

We determine a rest wavelength of 34.815 ± 0.004pm. The line widths are 
consistent with the nominal instrumental resolution except at P3 (ionization 
bar), where the FWHM is 88±3 km/s, less than the expected value of 1 1 4 km/s. 

Because the source is imaged on the individual field mirrors in the CGS (Erickson 
et al. 1985), a possible explanation is that the bulk of the [Sill] line 
intensity at P3 originates in an area significantly smaller than the beam. The 
apparent continuum at P6 has a definite slope not seen elsewhere. This spectrum 
was fitted using the slope of the continuum as a fifth free parameter. The FWHM 
and the wavelength of the fitted line profile are consistent with those found at 
the other positions. The line intensities are plotted versus position in Fig. 3. 

III. INTERPRETATION AND DISCUSSION 

Here we present some simple calculations which constrain the physical 
conditions in the [SiII](34.8pm) emitting region and then compare the observed 
intensities with the results of more detailed calculations to estimate the Si + 
abundance. 

A. The Photodissociation Region ( PDR) 

Observational evidence and theoretical considerations indicate that, with 
the exception of the localized ( —1 ' ) emission north of BN-KL, all of the extended 
[SiII](34.8pm) emission arises from the PDR in Orion, and not from the molecular 
cloud as a whole or from shocked gas. Gbservationally, the [Sill] emission 
spatially correlates with the extended [OI](63» 145pm) and [CII](158pm) emission 
(Ellis and Werner 1985), which have already been successfully modelled as arising 
in the PDR (Tielens and Hollenbach 1985b). Theoretically, one would not expect Si + 
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to exist in the dark molecular interior of OMC-1 , which is shielded from the 
photons which ionize Si. Alternatively, the silicon is predominantly in higher 
ionization states in the H1I region surrounding 0 l C Ori (Rubin, private 
communication). In PDRs, however, nearly all of the gas phase silicon is 
expected to be in the form of Si + . 

Face-on PDRs such as Orion are expected to have insufficient column densities 
to produce optically thick [Sill] transitions. If all the PDR gas-phase silicon 
lies in the ground state of Si + , the maximum line-center optical depth of [Sill] 
is T max * N 22 6 /Av, where Av is the line-of-sight velocity dispersion in units of 
km/s, 6 is the gas phase silicon abundance relative to its solar abundance x 

© 

(Table 1), and N 22 is the hydrogen column density in units of 10 22 cm -2 

(Hollenbach and McKee 1979). The grain attenuation of ultraviolet photons limits 

the extent of the warm gas in PDRs to A y < 3 or n 22 <. 0.6, 6 is probably £0.1 

(the value in diffuse gas), and velocity dispersions in PDRs are typically Av - 

1-5 (Tielens and Hollenbach 1985a,b), so that t <0.06. 

max 

Fig. 4 plots the intensity I of the [SiII](34.8pm) transition as a function of 
the column density N(Si + ) for various values of the hydrogen density n 0 and 
temperature T, assuming Av = 2. The linear dependence of I on N(Si + ) for values 
<, 3x10 * 7 cm~2 marks the optically thin regime. Using Fig. 4 and assuming N 22 < 

1, the observed value of I » 6 . 1 x 10*3 erg cm ~2 a -1 sr -1 at the Trapezium 
immediately shows that n 0 >, 10** cm - ^, T >, 100K, and 6 >, 0.03. 

Tielens and Hollenbach (1985b) have constructed a detailed model of the PDR 
behind the Trapezium, tailored to fit the observed far-infrared continuum and the 
[OI](63um, 145pm), [CII](158pm) , and [CI](609pm) integrated line intensities. 
Assuming 6 = 0.022 (the depletion determined from UV absorption measurements of 
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the heavily depleted cloud toward £ Oph; Morton 1974), they predicted an 
intensity of 1.8x10~3 ergs cnT^ sr~^ for the [SiII](34.8pm) line towards the 
Trapezium. This value is a factor of 3.4 smaller than we observe. The intensity 
of the [Sill] line scales with 6 in their model, so our observation implies 6 = 
0.075 or a silicon abundance of 2.6x10"^. These models are characterized by 
densities n 0 ~10 5 cnT^, temperatures T-500K, and column densities of warm hydrogen 
N 2 2 = 0.5 consistent with Fig. 4. 

The [Sill] intensity is sensitive to the collisional rates which are 
uncertain by ±50$ (cf Launay and Roueff 1977b with Dalgarno and McCray 1972). 

The Orion models of Tielens and Hollenbach (1985b) produce adequate fits (factor 
of 2) to the observations with peak PDR temperatures ranging from 250 K to 1000 K 
and n 0 -1~3x10 5 cm~3. including the absolute calibration uncertainty, we estimate 
that the silicon abundance is determined to within a factor of 2-3 and conclude 
that silicon is certainly no more depleted and probably somewhat less depleted in 
the denser (n o -10 s cm~3) Orion than in the less dense (n o <10 3 cm~3) £ oph. This 
is the first determination of silicon depletion in the outer layers of a dense 
molecular cloud. The trend toward increasing silicon depletion with increasing 
cloud density seen in diffuse clouds by UV absorption measurements (Snow and 
Joseph 1985) does not seem to extend into this dense PDR. 

B. The Shocked Region 

The emission peak north of BN could be produced by density or temperature 
fluctuations in the PDR gas. However, the spatial coincidence of this [Sill] 
peak with the peaks in the shock emission from H 2 , CO, [01], and OH is strong 
evidence that the enhanced [Sill] emission arises from shocked gas. The 
coincidence with the [01] shocked gas (Werner et al. 1984) is especially 
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convincing, since the [01] emission demonstrates the presence of a significant 
component of shocked atomic gas of similar excitation. Comparing the intensity 
at position P6 with that at the other map positions in Fig. 3, we estimate that 
the shock-produced emission is -4-7x1 0“^ ergs cnT^ s~^ sr -1 . 

Hollenbach (1982) and Chernoff et al. (1982) model the shocked gas in BN-KL 
as consisting of two components: a shocked outflow or wind with shock velocities 
v -100 km/s sufficient to dissociate and ionize the preshock gas and create a "J" 
shock, and a molecular "magnetic” or "C" shock caused by previously shocked 
supersonic gas sweeping up ambient molecular gas (see Dralne 1980 or McKee et 
al. 1984). Hollenbach (1985) shows how the [OI](63pm) luminosity from the wind 
shock might be used to measure the mass loss rate from the protostar, and Werner 
et al. (1984) model the [01] shock emission in BN-KL as a wind shock, obtaining M 
- 3x1 0”3 M©/yr. They argue that the [01] emission is too intense to be explained 
by the molecular shock, which has been invoked by Draine and Roberge (1982) and 
Chernoff, Hollenbach, and McKee (1982) to model the H 2 , CO, and OH emission. 
However, the C shock model of Draine and Roberge comes within a factor of a few 
of matching the observed [01] intensity if a large fraction of the oxygen is 
assumed to be in the atomic form in the preshock molecular gas. Thus, the 
analysis of the [OI](63pm) emission does not conclusively establish the existence 
of a wind J shock. 

The observation of shocked [Sill] provides additional evidence that the [01] 
and [Sill] emission arises from a wind J shock and not the molecular C shock. In 
the C shock, the hydrogen is mostly molecular and any pre-existing Si + is quickly 
converted to SiH + and SiH 2 + by reactions of "hot" Si + and SiH + ions with H 2 . On 
the other hand, the v - 100 km/s J shock provides a natural explanation for the 
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[Sill] emission. This shock dissociates and ionizes any pre-existing molecular 

gas phase silicon (cf Shull and McKee 1979). In addition, the shock sputters 

solid phase silicon from grains Into the gas (cf Draine and Salpeter 1979, Seab 

and Shull 1983). The silicon is ionized to Si ++ and Si +++ in the hot (T ~ 10 s K) 

post shock gas, but as the gas cools to T < 10 4 K downstream from the shock 

front, the Si +++ and Si ++ recombine to Si + . The recombination of Si + is offset 

by the photoionization caused by the FUV photons from the hot postshock gas. The 

hydrogen column density of warm (T >. 100 K) FUV-ionized gas can be as high as 

10 21 cm -2 if the preshock density is n >,10 5 cm* 3 (see Hollenbach 1980). Thus, 

ps * 

if nearly solar abundances of Si + are assumed and n > 10 s cm* 3 . a column N(Si + ) 

ps 

= 3*10* 6 cm~ 2 can be obtained in the J shock model, consistent with the 
observational results (cf Fig. 4). 

These qualitative conclusions are demonstrated in a more quantitative 
fashion by the detailed results of numerical J shock calculations (McKee, 
Chernoff, and Hollenbach 1984). These results include the destruction of 
preshock silicate grains and predict [SiII](34.8pm) intensities of 2.4x1 o _i< 
ergs cm -2 s -1 sr -1 for n pg = io 1 * cm* 3 and v g = 100 km/s, and 7.6x10* 3 
ergs cm -2 s" 1 sr _1 for n pgj = io 6 cm* 3 and v g = 100 km/s. Chernoff et al^ (1982) 
estimated n pg ~ io 5 cm* 3 in the wind shock around IRc2; therefore the wind shock 
model naturally provides an explanation for the silicon emission from the shocked 
gas. The observed intensity of shocked [SiII](34.8pm) is sufficiently high to 
suggest that either preshock silicate grains hqve been destroyed in the wind 
shock or that the preshock wind carries a near solar abundance of gas phase 
silicon. 

In conclusion, the [SiII](34.8ym) line is a valuable diagnostic of the 
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physical conditions in shocked atomic gas and the grain destruction processes 
which may occur therein. It is also an important new probe of the gas phase 
silicon depletion in dense (-10 s cm~3) interstellar clouds which cannot be 
observed with conventional UV techniques. 
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Table 1 


Strong Cooling Lines in Warm Atomic Gas 


Species \ hv/k I.P. a x b N„(t=* 1) C n b 

© H ^ p cr « 

(ym) ( °K) (eV) (cm ) (cm 


CII 

158 

92 

11.3 

3.7(-4) 

1 .2(21 ) 

2.8(3) 

01 

63 

228 

13.6 

6.7(-4) 

1.0(21) 

4.7(5) 

Sill 

35 

in 3 

8.2 

3.5(-5) 

4.6(23) 

3.4(5) 


cL k 

Ionization potential of the neutral species. D Solar abundances relative to 
hydrogen (Allen 1973). c Column density of hydrogen required for optical depth 
unity at line center, assuming a line-of-sight velocity dispersion of Av = 1 km/s 
and a solar abundance of the species in its ground state. Critical atomic 
hydrogen density; collisional rates for CII from Launay and Roueff (1977b), for 
Oil from Launay and Roueff (1977a), and for Sill from Dalgarno and McCray (1972). 
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Figure Captions 


Figure 1. The beam locations at which [SiII](34.8pm) was observed. 

Starting in the lower right, the three crosses mark the positions of © 2 0ri A, 
0*Ori C, and the far-infrared continuum peak of the Kleinmann-Low nebula. 

Figure 2. The [SiII](34.8pm) line and neighboring continuum (a) at the 
ionization bar (P3)> (b) at Q*Ori C (P5), and (c) NW of the Kleinmann-Low nebula 
(P6). 

Figure 3. The [SiII](34.8pm) line intensity as a function of position 
across Orion-KL. The error bars are one standard deviation of the mean; they 
represent formal statistical errors only. 

Figure 4. A theoretical prediction of the [SiII](34.8pm) line intensity 
versus the column density of Si + for different combinations of density and 
temperature (n 0 , T). n 0 » ® corresponds to n 0 » n » 3.4x10 s cm - ^; T = » 
corresponds to T >> 413 K. cr 
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[Sill] emission from both the extended photodissociation region surrounding 
3 1 0ri C and from the shocked gas NW of BN-KL. The inferred gas-phase silicon 
element al abundance relative to hydrogen in the dense (~10 5 cm -3 ) primarily 
neutral photodissociation region is approximately 2.6 x 10 -6 , a factor of 
3.075 times the solar value and 3.4 times greater than the abundance in the 
noderate density (<10 3 cm -3 ) cloud toward c, Oph. The silicon abundance in 
-he shocked gas is approximately solar, indicating that any pre-existing grains 
nave been destroyed in the shock wave or that the preshock gas carries a near 
solar abundance of gas phase silicon. The shock-excited [Sill] (34.8 pm) 
emission may arise from shocked wind material in the outflow around IRc2, with 
•find velocities >100 km/s. 
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